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A SIMPLIFIED METHOD FOR CALCULATING TEMPERATURE 
TIME HISTORIES IN CRYOGENIC WIND TUNNELS 

Robert L. Stallings, Jr. 

Langley Research Center 

Mi 1 ton Lamb 

Langley Research Center 
SUMMARY 

A simplified method for calculating average temperature time histories of 
the test gas and tunnel walls for cryogenic wind tunnels has been developed. 
Results from the method are in general agreement with limited preliminary 
experimental measurements obtained in the Langley 13.5-Inch P^’lot Cryogenic 
Wind Tunnel. 


INTRODUCTION 

The advantages of operating a wind tunnel at cryogenic temperatures in 
order to obtain high test Reynolds number has been well documented in recent 
years. (See, for example, refs. 1 and 2.) The major advantages of operating 
at an extremely low temperature are reduced dynamic pressure for a given 
Reynolds number and, for fan-driven tunnels, greatly reduced drive-power 
requirements as compared to conventional wind tunnels. The need for improved 
Reynolds number capability, the advantages of a cryogenic tunnel with regard 
to. capital cost, and the current national need to reduce energy consumption 
(ref. 3) were factors which contributed to the United States Government's 
decision to construct a large transonic tunnel at the Langley Research Center 
capable of operating at cryogenic temperatures. This facility, to be known as 
the National Transonic Facility, will have a 2.5 by 2.5 meter test section and 
will have a Reynolds number capability based on typical model chord of 
approximately 120 x 10^ when operating in the cryogenic mode at its maximum 
operating pressure of 8.8 atmospheres. 
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Some early preliminary studies conducted at the Langley Research Center 
related to the feasibility of operating at cryogenic temperatures were 
concerned with the time required for tunnel "cool down" for cryogenic operation, 
the amount of coolant required to maintain the cryogenic temperatures during 
a cryoqenic run, and the time required for tunnel "warm up" after completion of 
a cryogenic run. In order to estimate these parameters for various insulation 
thicknesses, structural weights, and test gas properties, a simplified method 
for calculating average conditions for the entire tunnel system was developed. 
This method is described in this report and some results are compared with 
preliminary measurements obtained in the Lanqley 13.5'Inch Pilot Cryoqenic 
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SYMBOLS 

2 

surface area of tunnel interior wall, m 

2 

test section cross-sectional area, m 

width of air gap, m 

specific heat of test gas, ^ 

average heat transfer coefficient for tunnel interior wall, 
J 

2 

m - sec - K 


average heat transfer coefficient for natural convection over 

tunnel exterior wall, 

nr - sec - K 

thermal conductivity of air, - — — 7 

n — sec ~ i\ 


coolant mass flow rate. 


_k_S_ 

sec 


test section Mach number 


stagnation pressure. 
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q 


AMB 


w,Bw 


1 1 1 rs 
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test section static pressure, ~ 

ITT 

heat flux, ^ 
time, sec 

atmospheric temperature, K 
saturation temperature of coolant, K 
test gas temperature, K 

tunnel interior wall surface temperature, K 
temperature of external wall surface, K 
wall temperatures across air gap, K 


w mass of test gas, kg 

\ Latent heat of vaporizaiion of coolant, - j^ -— 

r power coefficient ( 0.2 for 0.1! My^ll.O) 

ANALYSIS 

The simple heat balance used in this analysis is based on average 
conditions throughout the tunnel system. For example, heat convected from the 
system is based on an average heat transfer coefficient, conduction through the 
tunnel walls is based on an average wall thickness, and the test gas 
temperature is calculated as an average temperature. With this in mind the 
heat balance equation becomes 


q test gas - q N 0 " q conv. 


xomp , 


( 1 ) 


where q. . is the heat stored in the test gas, q M is the cooling 

test gas in ^ 

available from the liquid nitrogen, q conv is the heat exchanged with the 
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tunnel w,» 1 I , and q c&njn is the energy flux transferred to the test gas by 
the drive fnn. Expressing the heat flux parameters in terms of the average 
spacial temperatures and thermal properties (the expression for q comn is 
from ref. 2), equation (1) becomes: 

"l - * ("» + c p (I T - T s ) -hA(T w - T t ) 

- 14.05 P TS M 3 ts VT t A ts 1 

where T t is the average test gas temperature at time t. Ihe variation of 
Ty with time was obtained by solving equation (2) using a finite difference 
scheme resulting in the following equation: 


' T r) + 

Tn 
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Therefore, if conditions at some initial time t are known the tunnel test gas 


temperature after a time increment t can be calculated. After the tunnel 


test gas temperature is calculated at time t + ‘t the change in tunnel wall 
temperature for this time increment can be calculated. For the present method, 
the wall temperatures were calculated from the program of reference 4 modified 
to include equation (3). 

The finite difference solution of reference 4 applies to the case of a 
thermally thick wall with various modes of heat transfer, including convection, 
occurring at the exterior surface; however, the interior surface or back wall 
is assumed to be insulated. When applied to the tunnel wall case, the exterior 
surface as defined in reference a is actually the tunnel wall interior surface 
whereas the back wall as defined in reference 4 would be the exterior surface • 
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of the tunnel wall woich is exposed to atmospheric or room condition. In tnis 
discussion terminology applicable to the tunnel wall case is used. Since an 

exchange of heat from the exterior wall and its surrounding environment will 
definitely occur, it is necessary that it be accounted for in the calculation. 
The desired noninsulated exterior wall boundary conditions were artifically 
ri btained by making the proper choice of variables and modes of heat transfer 
available as options in the program that would be representati ve of the actual 
transfer of heat. 

Two noninsulated conditions were considered for the present method. First 
and simplest of the two is the case where the exterior wall temperature is held 
constant at ambient room temperature. This condition was simulated in the 
program bv locating a fictitious block having a large thermal mass adjacent to 
the exterior wall surface. The conduction o f heat to his fictitious block, 
maintains the exterior wall surface temperature at an essentially constant 
value and therefore correctly simulates the flow of heat through the tunnel 
wall for this hypothetical constant temperature boundary condition. 

The second noninsulated condition considered for the exterior wall is 
the case of free convection to a standard atmospheric environment. The flow 
of the heat for this condition is governed by the equation 

q = h F A(T w,Bw * T AMB ) i4) 

where hp is the coefficient of heat transfer for free convection conditions, 

T w Bw is the temperature of the exterior wall, and T^ is the ambient 

atmospheric temperature. Equation (4) is very similar in form to the following 
equation (eq. (7), ref. 4): 
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for heat conduction across an air gap. The symbol k . is the conductivity 

31 r 

of air, T , and T 0 are the wall temperatures of the blocks on each side of 

W f I W f t 

the gap, and b is the width of the gap or the distance between the two blocks 
1 and 2. If k . / b is set equal to h F , T , set equal to T R , and T 0 

3li I W y I W 9 BW W^l 

set equal to t aMB ’ then ec l uat10ns (4) and (5) are equal. If b is assumed 
to have unity value, then the value of k . will equal the value of h c 
Therefore, free convection at the exterior wall can be artificial l^imulated in 
the program by assuming an air gap separates the exterior wall from a fictitious 
block having a very large density and specific heat. The initial temperature of 
the fictitious block is set equal to the ambient temperature and its mass and 
specific heat are fixed at sufficiently large values such that it remains at 
ambient temperature for the maximum time span of interest. The tunnel wall 
temperatures computed from the program would then represent the proper values 
for the case where the exterior wall is exchanging heat by free convection to 
an atmospheric environment. 


RESULTS 

In order to determine if the present method gave realistic estimates of 
temperature-time histories for the operation of a cryogenic tunnel, comparisons 
were made between calculations and experimental measurements of the test gas 
temperature obtained in the Langley 13.5-Inch Pilot Cryogenic Wind Tunnel 
(ref. 5). Tunnel speci ficiations that were used as inputs to the program are 
shown in Table 1. Calculations were made for the pre-run cool-down period and 
the post-run warm-up period. During these periods, the exterior wall 
temperature was assumed to remain constant at 294 K and was incorporated into 
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the program using the first noninsulated method previously discussed. 

Typical calculated cool-down temperature time histories are shown in 
figure 1 for M.^ = 0.4 and P^. = 1 atmosphere. Temperatures are shown for 
the test gas, the interior wall, and the exterior wall ( 294 K). The 
calculations indicate that the test gas temperature and interior wall 
temperature decrease approximately linearly with time and reach the desired 
78 K operating temperature in approximately 30 minutes. Also shown in figure 1 
is a measured data point for the facility from reference 5 for similar test 
co r ditions. This measurement indicates the time required to cool the facility 
from aoproximately 294 K to 106 K is approximately 30 minutes as compared to 
apDroximately 24 minutes predicted by the Dresent method. 

Calculated temperature distributions through an average tunnel wall 
section are shown in figure 2 for several values of t during the ccol-down 
period. Due to the large thermal -conductiv ity of aluminum, the temperature is 
essentially constant from the interior wall to the aluminum-urethane interface. 
Large temperature gradients occur through the urethane with the largest gradients 
occurring at the aluminum-urethane interface. 

Shown in figure 3 is the effect of test section Mach number on calculated 
cool-down temperature-time histories for the test gas at a constant tunnel 
total pressure equal to 1 atmosphere. At the higher test-section Mach numbers, 
the results indicate a rapid increase in time required for cool-down with 
increasing Mach number. Although cool-down time measurements have not been 
obtained in the pilot tunnel over a wide range of Mach number, the trend of 
increasing cool-down time with increasing Mach number has been observed in 
the facility. This increase in ‘"’’me is a direct result of the increasing 
energy flux transferred by the fan to the stream with increasing Mach number. 
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As shown in equation (3), this energy varies as H 

Shown in figure 4 are calculated temperature-time histories for the test 
gas in the Pilot Cryogenic Tunnel during warm-up periods for stagnation 
pressures of 1 and 2 atmospheres and a test section Mach number of 0.6. During 
the warm-up perio' 4 which occurs after completion of a cryogenic run, the flow 
of liquid nitrogen is terminated and the flow of heat through the tunne 1 walls 
and from the drive-fan increases the test gas temperature as indicated in 
figure 4. The calculated warm-up time for the case of 2 atmospheres tunnel 
pressure can be compared with measurements from reference 5 that are shown in 
figure 4 for similar conditions. The measured warm-up time from 1:5 K to 300 K 
is aDproximately 53 minutes as compared to approximately 60 minute;' calculated 
f^cm the present method. The calculations indicate that the warm-up time varies 
inversely proportional to the tunnel stagnation pressure. 

The general agreement shown between calculations and measurements in 
f igures 1 and 4 indicates that the present method does give realistic trends 
for the variation of average temperatures in a cryogenic wind tunnel. It 
should be recognized, however, that the results are for average conditions and 
can vary significantly from local conditions throughout the tunnel. 

CONCLUDING REMARKS 

A simplified method for calculating average temperature -time histories of 
the test gas and tunnel walls for cryogenic wind tunnels has been developed. 
Results from the method are in general agreement with limited preliminary 
experimental measurements obtained in the Langley 13.5-Inch Pilot Cryogenic 
Wind Tunnel . 
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Table 1. - PROGRAM INPUTS FOR 13. 5- INCH PILOT CRYOGENIC WIND TUNNEL 

1) Structural Material, Aluminum 

(a) Mass = 3175kg 

2 

(b) Surface Area = 55.74 m 

/ \ « t u * l _ _ Mass = 2.11 cm 

(c) Average Thickness = — 

Surface Area x Density 

2) Insulation (External), Urethane 
Average Thickness = 12.7 cm 

3) Volume = 17 m^ 

4) Coolant, Liquid N, 

(a) Z 78 Y. Flow rate = 1.02 ^ 

(b) T^ = 78 K: Flow rate regulated to maintain this temDerature. 
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Figure 1.- Typical temperature-time histories during cool-down of Langley 
13.5-Inch Pilot Cryogenic Wind Tunnel calculated from present 

method, Mys = o o 4, Pt = K0 atm„ 
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Figure 4 C - Effect of tunnel pressure on calcu 
5 histories during warm-up of 13.5- 




